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We study the onset of a classical order parameter after a second-order phase transition in quantum 
field theory. We consider a quantum scalar field theory in which the system-field (long-wavelength 
modes), interacts with its environment, represented both by a set of scalar fields and by its own 
short- wavelength modes. We compute the decoherence times for the system-field modes and compare 
them with the other time scales of the model. We analyze different couplings between the system 
and the environment for slow quenches. Within our approximations decoherence is in general a 
short time event. 

PACS numbers: 03.70.+k, 05.70.Fh, 03.65.Yz 

I. INTRODUCTION 

The emergence of classical physics from quantum behaviour is important for several physical phenomena in the 
early Universe. This is beyond the fundamental requirement that only after the Planck time can the metric of the 
Universe be assumed to be classical. For example, 

• The inflationary era is assumed to have been induced by scalar inflaton fields, with simple potentials [1]. Such 
fields are typically assumed to have classical behaviour, although in principle a full quantum description should 
be used. 

• The origin of large scale structure in the Universe can be traced back to quantum fluctuations that, after crossing 
the horizon, were frozen and became classical, stochastic, inhomogeneities [2]. 

• It is generally assumed that several phase transitions have occurred during the expansion of the Universe [3] . As 
in the case for the inflaton fields, the (scalar) order parameter fields that describe these transitions are described 
classically. However, the description of early universe phase transitions from first principles is intrinsically 
quantum mechanical [4]. 

• As a specific apphcation [5] of the previous point, the very notion of topological defects (e.g. strings and 
monopoles) that characterize the domain structure after a finite-time transition, and whose presence has con- 
sequences for the early universe, is based on this assumption of classical behaviour for the order parameter [6], 
as it distributes itself between the several degenerate ground states of the ordered system. 

In the present paper we are concerned with the third point above, the quantum to classical transition of the order 
parameters in second order phase transitions. Any approach must take into account both the quantum nature of the 
order parameter and the non-equilibrium aspects of the process. The problem of the quantum to classical transition 
in the context of inflationary models was first addressed by Guth and Pi [7]. In that work, the authors used an 
inverted harmonic oscillator as a toy model to describe the early time evolution of the inflaton, starting from a 
Gaussian quantum state centered on the maximum of the potential. They subsequently showed that, according to 
Schrodinger's equation, the initial wave packet maintains its Gaussian shape (due to the linearity of the model). 
Since the wave function is Gaussian, the Wigner function is positive for all times. Moreover, it peaks on the classical 
trajectories in phase space as the wave function spreads. The Wigner function can then be interpreted as a classical 
probability distribution for coordinates and momenta, showing sharp classical correlations at long times. In other 
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words, the initial Gaussian state becomes highly squeezed and indistinguishable from a classical stochastic process. 
In this sense, one recovers a classical evolution of the inflaton rolling down the hill. 

A similar approach has been used by many authors to describe the appearance of classical inhomogeneities from 
quantum fluctuations in the inflationary era [8]. Indeed, the Fourier modes of a massless free fleld in an expanding 
universe satisfy the linear equation 

cl>'i + (fc^ - ^)</.fe = 0. (1) 

For sufSciently long- wavelengths (fc^ ^ a" /a), this equation describes an unstable oscillator. If one considers an 
initial Gaussian wave function, it will remain Gaussian for all times, and it will spread with time. As with the toy 
model of Guth and Pi, one can show that classical correlations do appear, and that the Wigner function can again be 
interpreted as a classical probability distribution in phase space. [It is interesting to note that a similar mechanism can 
be invoked to explain the origin of a classical, cosmological magnetic field from amplification of quantum fluctuations] . 

However, classical correlations are only one aspect of classical behaviour. It was subsequently recognized that, in 
order to have a complete classical limit, the role of the environment is crucial, since its interaction with the system 
distinguishes the field basis as the pointer basis [9]. [We are reminded that, even for the fundamental problem of the 
space-time metric becoming classical, simple arguments based on minisuperspace models suggest that the classical 
treatment is only correct because of the interaction of the metric with other quantum degrees of freedom [10].] 

While these linear instabilities cited above characterise free fields, the approach fails when interactions are taken into 
account. Indeed, as shown again in simple quantum mechanical models (e.g. the anharmonic inverted oscillator), an 
initially Gaussian wave function becomes non-Gaussian when evolved numerically with the Schrodingcr equation. The 
Wigner function now develops negative parts, and its interpretation as a classical probability breaks down [11]. One 
can always force the Gaussianity of the wave function by using a Gaussian variational wave function as an approximate 
solution of the Schrodingcr equation, but this approximation deviates significantly from the exact solution as the wave 
function probes the non-linearities of the potential [11,12]. 

When interactions are taken into account, classical behaviour is recovered only for "open systems", in which the 
relevant degrees of freedom interact with their environment. When this interaction produces both a diagonalization 
of the reduced density matrix and a positive Wigner function, the quantum to classical transition is completed [13]. 

Going from quantum mechanical toy models to quantum field theory is, of course, extremely difiicult. For this 
reason, several authors [14] have considered different approximations in quantum mechanics, and compared them to 
the exact results. If successful in quantum mechanics, they could be implemented with greater confidence in field 
theory. As already mentioned, when this procedure is applied in the context of closed systems^ the conclusion is 
that, in general, mean-field approximations do not reproduce the evolution of the system at late time [11,12,14]. 
Therefore, in principle, there is not reason to believe they will do so in field theory. In spite of this, computational 
necessity has lead many authors to perform calculations in closed field theories with Hartrce, mean field, or 1/A'^ 
approximations, since they are non-perturbative, well-defined and suitable for numerical calculations [15-18]. In such 
calculations, classical correlations do appear in some field theory models [16,17]. However, since such decoherence (in 
a time-averaged sense) takes place at long times after the transition has been achieved initially, when the mean field 
approximation has broken down, this may be an artifact of the Gaussian- like approximations [11]. 

In a previous paper [11] we considered similar arguments for quantum mechanics, but for open systems. When 
an anharmonic inverted oscillator is coupled to a high temperature environment, it becomes classical very quickly, 
even before the wave function probes the non-linearities of the potential. Being an early time event, the quantum to 
classical transition can now be studied perturbatively. In general, rccohcrence effects are not expected [19]. Taking 
these facts into account, we have extended the approach to field theory models [20]. In field theory, one is usually 
interested in the long- wavelengths of the order parameter. Moreover, the early universe is replete with fields of all sorts 
which comprise a rich environment. For this reason, we considered a model in which the order parameter interacts 
with a large number of environmental fields, including its own short- wavelengths. Assuming weak coupling and high 
critical temperature, we have shown that decoherence is a short time event, shorter than the spinodal time tgp, which 
is essentially that time by which the order parameter field has sampled the degenerate ground states. As a result, 
perturbative calculations arc justified [20]. Subsequent dynamics can be described by a stochastic Langevin equation, 
the details of which are only known for early times. 

Our approach in Ref. [20] has some connections with well-established classical behaviour of thermal scalar field 
theory [21] at high temperature. In many articles it has been shown that, at high temperatures, the behaviour of 
long-wavelength modes is determined by classical statistical field theory. The effective classical theory is obtained 
after integrating out the hard modes with k > gT. Although similar, this approach has some important differences 
from ours: the "classical behaviour" in this soft thermal mode analysis is defined through the coincidence of the 
quantum and the statistical correlation functions. In particular, thermal equilibrium is assumed to hold at all times. 
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Finally, the cutoff that divides system and environment depends on the temperature, which is externally fixed. In our 
approach, the quantum to classical transition is defined by the diagonalization of the reduced density matrix, which 
is not assumed to be thermal. In phase transitions the separation between long and short-wavelengths is determined 
by their stability, which depends on the parameters of the potential. 

In Ref. [20] we only considered the case of an instantaneous quench, and bi-quadratic coupling between the system 
and the environment. In Ref. [22] we began to extend those results to the case of a slow quench. In this paper 
we provide the details of that analysis and extend it to other couplings between system and environment. The 
consideration of slow quenches is very important since the Kibble-Zurek mechanism predicts the relation between the 
subsequent domain structure and the quench time [23-26] (by indirectly counting defects). 

The paper is organized as follows. In Section 2 we introduce our models. These are theories containing a real 
system field (f), which undergoes a transition, coupled to other scalar fields Xa (a = 1)2, which constitute the 

external part of the environment. Gauge fields bring their own specific difficulties and we shall not discuss them here. 
We compute the influence functional by integrating out the environmental fields for different couplings. Section 3 is 
dedicated to reviewing the evaluation of the master equation and the diffusion coefficients which are relevant in order 
to study decoherence. In Section 4 we evaluate upper bounds on the decoherence time for slow quenches. As we will 
see, provided quenches are not to slow, decoherence takes place before the field samples the minimum of the potential, 
i.e. decoherence time is typically shorter than the spinodal time. However, since all relevant timescales depend only 
logarithmically upon the parameters of the theory, it is necessary to keep track of 0(1) prefactors carefully, something 
that we rather took for granted in Ref. [20] and in Ref. [22]. We will also show that the bi-quadratic coupling is 
the most relevant for the quantum to classical transition. Section 5 contains the conclusions of our work. Two short 
appendices fill in some of the detail. 



II. THE INFLUENCE ACTION FOR AN EXTERNAL ENVIRONMENT 



For the infinite degree of freedom quantum field (f> undergoing a continuous transition, the field ordering after 
the transition begins is due to the growth in amplitude of its unstable long-wavelength modes. For these modes 
the environment consists of the short-wavelength modes of the field, together with all the other fields with which 4> 
inevitably interacts [27-29] in the absence of selection rules. The inclusion of explicit environment fields is both a 
reflection of the fact that a scalar field in isolation is physically unrealistic, as well as providing us with a systematic 
approximation scheme. 

The (/i-field describes the scalar order parameter, whose Z2 symmetry is broken by a double-well potential. Specif- 
ically, we take the simplest classical action with scalar and environmental fields Xa 

'S' [</>, X] = 'S'syst [</»] + -Senv [x] + 'S'int [Xa , 0] , (2) 

where (with fj,"^, w? > ) 



SeM = jzj {^^^Xa^'^Xa - ^m^X^j . (4) 



The most important interactions will turn out to be of the biquadratic form 

N 

5int[Xa,<A] = S^4<j),x] =-IljJ d'x cl,\x)xl{x). (5) 

Even if there were no external x fields with a quadratic interaction of kind (5), the interaction between short and 
long- wavelength modes of the ^field can be recast, in part, in this form (see later), showing that such a term is 
obligatory. The generalization to a complex field (p is straightforward, and has been considered elsewhere [22] . 
Later, we shall consider additional interactions to the biquadratic interaction of a bilinear form 

N I 

5biiin[</', x] = -Y.^-f d'^x 4>\x)x.{xl (6) 

a=l 
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in which the environment couples hnearly. 

We could also have included fermionic Yukawa interactions but these provide a smaller contribution to the diffusion 
constant because of Fermi statistics (which gives a diffusion coefficient relatively 0(/x^/r^)). However, because the 
biquadratic self- interaction is the overwhelming term in the diffusion coefficient, in general we are not interested in 
effects that are overshadowed by it, and we shall not consider Yukawa interactions here. 

Our exception to this rule of ignoring small contributions is in the the inclusion of linear couplings of the form 

^ It 2 r 

S^^[4>, Xa] = - E ^ / '^'^ <^(^)Xa(a;), (7) 

a=l •' 

in which the order parameter field couples linearly to a linear environment. Much if not most of the work on 

decoherence has been for linear coupling to the environment of this type (e.g. see [30,31]). While it is sensible 
in quantum mechanics, in the context of quantum field theory such linear couplings signal an inappropriate field 
diagonalisation, although they are exactly solvable in some circumstances [31]. 

Although the system field (p can never avoid the decohering environment of its own short-wavelength modes, to 
demonstrate the effect of an environment we first consider the case in which the environment is taken to be composed 
only of the fields Xa- We are helped in this by the fact that environments have a cumulative effect on the onset of 
classical behaviour. That is, the inclusion of a further component of the environment reduces the time i/j it takes for 
the system to behave classically (for the definition of see later). Thus it makes sense to include the environment 
one part after another, since we can derive an upper bound on to at each step. The short-wavelength modes of the 
field will be considered last. 

To keep our calculations tractable, we need a significant part of the environment to have a strong impact upon the 
system- field, but not vice- versa, from which we can bound t/j. The simplest way to implement this is to take a large 
number iV ^ 1 of scalar Xa fields with comparable masses ma — At weakly coupled to the (j), with A, Qa, 9^,, s'J. 1- 
Thus, at any step, there are N weakly coupled environmental fields influencing the system field, but only one weakly 
self-coupled system field to back- react upon the explicit environment. 

We first consider the case in which the fields Xa interact through the biquadratic interaction (5) alone. For one-loop 
consistency at second order in our calculation of the diffusion coefficient (that enforces classicality) it is sufficient, 
at order of magnitude level, to take identical ga, = g/\fN. Further, at the same order of magnitude level, we take 
g — X. This is very different from the more usual large-N 0{N + l)-invariant theory with one ^-field and N Xa fields, 
dominated by the 0(1/ N) (x^)^ interactions, that has been the standard way to proceed for a closed system. With 
our choice there are no direct x'^ interactions, and the indirect ones, mediated by ^ loops, are depressed by a factor 
g/VN. In this way the effect of the external environment qualitatively matches the effect of the internal environment 
provided by the short-wavelength modes of the (/)-field, but in a more calculable way. 

For small g the model has a continuous transition at a temperature T^. The environmental fields Xa reduce and, 
in order that T^/lJ.'^ = 24/(A + J29a)^ 1, we must take X + ^ga,<^l, whereby 1 > I/a/ZV » g. Further, with this 
choice the dominant hard loop contribution of the (^-field to the Xa thermal masses is 

Jml = 0{gT^/sfN) = 0{i^^/N) «; /i^. 

Similarly, the two-loop (setting sun) diagram which is the first to contribute to the discontinuity of the X-^eld 
propagator is of magnitude 

g^T^/N = 0{g„yNy^) « <5m|, 

in turn. That is, the effect of the thermal bath on the propagation of the environmental x-ficlds is ignorable. This 
was our intention; to construct an environment that reacted on the system field, but was not reacted upon by it to 
any significant extent. In particular, the infinite A'' limit does not exist. Depende nce on N is implicit through Tc as 
well as through the couplings, for initial temperatures Tq — 0{Tc). With rj = \J determining the position of 
the minima of the potential and the final value of the order parameter, this choice of coupling and environments gives 
the hierarchy of scales necessary for establishing a reliable approximation scheme, as has been shown in [20]. 

We shall assume that the initial states of the system and environment are both thermal, at a high temperature 
To > Tc. We then imagine a change in the global environment (e.g. expansion in the early universe) that can be 
characterised by a change in temperature from Tq to Tf < T^. That is, we do not attribute the transition to the effects 
of the environment-fields. 

Given our thermal initial conditions it is not the case that the full density matrix has and x fields uncorrelated 

initially, since it is the interactions between them that leads to the restoration of symmetry at high temperatures. 
Rather, on incorporating the hard thermal loop 'tadpole' diagrams of the x (and </>) fields in the (/> mass term leads 
to the effective action for ^ quasiparticles, 
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5.fstM = J d'x I ia^<^a''</, - ^m|(To)</.2 - A<^4| (8) 

1{Tq) cx (1 — Tq/Tc) for T « Tc- As a result, we can take an initial factorised density matrix at temperature 

To of the form p[Tq] = p,p[Tq] x p^[To], where /5^[T'o] is determined by the quadratic part of 5gyst[(/)] and /5^[To] by 
SenviXa]- That is, the many Xa fields have a large effect on (j), but the 0- field has negligible effect on the Xa- 

Provided the change in temperature is not too slow the exponential instabilities of the 0-field grow so fast that 
the field has populated the degenerate vacua well before the temperature has dropped significantly below T^. Since 
the temperature has no particular significance for the environment fields, for these! early times we can keep the 
temperature of the environment fixed at T^. « (our calculations are only at the level of orders of magnitude). 
Meanwhile, for simplicity the Xa masses are fixed at the common value m ~ /z. 

Our interests in the transitions in the early imiverse (domain structure, defect formation and, indirectly, structure 
formation) and the subsequent field evolution [6] are naturally couched in terms of fields, rather than particles. Since 
we need to be able to distinguish between different classical system-field configurations evolving after the transition, 
we will only be interested in the field-configuration basis for this reduced density matrix (in analogy with the usual 
quantum Brownian motion model, e.g. [32]). The resulting reduced density matrix [</)+, </)~, = ((/)+|/5r(t)|(/)~) 
describes the evolution of the system under the influence of the environment, and is defined by 

N 



= J ti'^Xa. p[(l>'^,XB.,(t> ,Xa.,t\, 



where p [</'''', X^ Xa' ^] = {4'~^xt\Pit)\4'~Xa) the full density matrix. For the reasons given above, the envi- 
ronment will have had the effect of making the system effectively classical once pi{t) is, approximately, diagonal. 
Quantum interference can then be ignored and we obtain a classical probability distribution from the diagonal part 
of Prit), or equivalently, by means of the reduced Wigner functional, which is positive definite after to [H]- For weak 
coupling there will be no 'recoherence' at later times in which the sense of classical probability will be lost [19]. 
The temporal evolution of Pi is given by 

Pr [(f>i ,(f)Y,t]= j d(j)'l j #r Jr [(t)f , (t)f ,t\cl)^, ^i" , to] Pr [(t>t > *o] , 

where Jr is the reduced evolution operator 



Jr[0+,^f-, to] = f ' [ ' Vr e^{^[^+]-^[^"l>T[<^+,r]- 

Jth+ Jd>7 



(9) 



The Feynman- Vernon [33] influence functional (IF) F[(/)^,(f> ] in (9) can be written in terms of influence action 
6A[(I)'^ and the coarse grained effective action (CGEA) A[0+,^~] by (see [28] for formal definitions) 

F[ct>+,r]=eKpi6A[cf>+,cf>-], (10) 



A[<p+,r] = s[<p+] - sir] + sA[<p+,r]- (n) 

Beginning from the initial distribution, peaked around </) = 0, we follow the evolution of the system under the 
influence of the environment fields. We will calculate the infiuence action to lowest non-trivial order (two vertices) for 
large N. Similarly to the biquadratic interaction, we assume weak coupling A ~ (/' ~ 5" <C 1, where we have defined 
g',g" by the order of magnitude relations g'^^ = g'/y/N and g'J^ = g"/y/N respectively. 

As we are considering weak coupling with the environment fields, we may expand the influence functional T [</>+, 0~] 
up to second non-trivial order in coupling strengths. The general form of the influence action is then [28,34] 

sA[<i>+,r] = {(5i„t[<^+,x^])o - (5int[r,Xa-])o} 

+ ^{(5ilt[</'+,X+])o- [(5i„t[0+,x^])o]'} 

-l{{Sir.t[^+,xt]Sint[r,X:]h - (5int[<^+,X^])o(^int[0-,Xa]>o} (12) 
+ l{{SL[r,Xa])0 - [(5i„t[r ,Xa]>0]'}. 
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For weak couplings it is relatively easy to compute the upper bound on the decoherence time due to the interactions 
Sint[(l>,x] of (5), (6) and (7). We shall find that, in our particular model, it is in general shorter than the spinodal 
time tsp, defined as the time for which 

{<t>'K -rf= 6/xVA . (13) 

[More exactly, t^-p should be defined as the time that the r.m.s of the field takes to reach the spinodal point, the point 
of inflection in the potential, but the difference is only logarithmically small, so we stay with (13).] 

It is for this reason that, with the qualifications below, we can use perturbation theory. In consequence, by the 
time that the field is ordered it can be taken to be classical. This has implications [22] for the formation of the defects 
that are a necessary byproduct of transitions. 



III. MASTER EQUATIONS AND THE DIFFUSION COEFFICIENTS 

In this Section we will obtain the evolution equation for the reduced density matrix (master equation), paying 

particular attention to the diffusion term, which is responsible for decoherence. We will closely follow the quantum 
Brownian motion (QBM) example [30,35], translated into ciuantum field theory [20,28]. 

This Section contains a gc;neralization of the; rc;sult already given in [20], adapted to the slow quenches that are 
physically relevant, in which m|(r) varies linearly in time. The formal deduction of the master equation is not 
altered by the duration of the quench. For the details we will follow Ref. [28] . The first step in the evaluation of the 
master equation is the c;alculation of the density matrix propagator Jp from Eq.(9). In order to solve the functional 
integration which defines the reduced propagator, we perform a saddle point approximation 

J,[(l)j,(l)-^,t\(j)tATM ~ expa[0+,(/)-i], (14) 
where 0^ is the solution of the equation of motion = with boundary conditions <^^(to) = (f>t 

Even then, it is very difficult to solve this equation analytically. We are helped by the observation that the ordering 
of the field is due to the growth of the long- wavelength unstable modes. Unstable long- wavelength modes start growing 
exponentially as soon as the quench is performed, whereas short- wavelength modes will oscillate. As a result, the 
field correlation function rapidly develops a peak (Bragg peak) at wavenumber fc = <C /x. Specifically [36], initially 
as hp' = 0{iJ,/ ^tTq), where is the quench rate. Assuming that a classical description can be justified post-hoc, a 
domain structure forms quickly with a characteristic domain size 0(k~^), determined from the position of this peak. 
[As an example, see the numerical results of [25] , where this classical behaviour has been assumed through the use of 
stochastic equations (see later)]. With this in mind, we adopt an approximation in which the system-field contains 
only one Fourier mode with k = Icq = 0(k~^), characteristic of the domain size. 

As "trial" classical solutions, we take 

0ci(f , s) = /(s, t) cos(fcoa;) cos(fcoj/) cos(fcoz), (15) 

where f{s,t) satisfies f{0,t) = (f>i and f{t,t) — 0f. Eq. (15) is an exemplary configuration that mimics domain 
formation. As intended, it represents domains of finite size k^^ . Although such a regular domain structure is idealised, 
numerical results (as in [25]) suggest that it is a reasonable first step. In practice, after this detailed motivation, we 
shall find that the decoherence time to is insensitive to wavelength for ko small. For the purpose of calculation, it 
is sufficient to set fco = 0, even though the physical situation requires a finite fco. An exact matching of fco to k is 
unnecessary. [We note that, in [20] and [22], we adopted a simpler form (j)c\{x, s) = f{s, t) cos{ko.x). Although showing 
a domain structure in one dimension is less physical, our conclusions are unchanged although details differ.] 
We write f{s, t) as 

f{s,t) = (l}iUi{s,t) + (j)tU2{s,t), (16) 

where, during the quench, Ui{s, t) are solutions of the mode equation with boundary conditions ui(0, t) — I, ui{t, t) = 
and 712(0, i) = 0, U2{t,t) = 1. In order to obtain the master equation we must compute the final time derivative of 
the propagator Jj.. After that, all the dependence on the initial field configurations (pf (coming from the classical 
solutions must be eliminated. The free propagator, defined as 

Jo[</.f+,<^r.*l<^^C>0] = V<l>+ I?r exp{i[5o(</.+) - 5o(r)]}; (17) 
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satisfies the general identities [28,35] (which is valid for instantaneous and for slow quenches) 



<l)f[u2{s,t) 



U2it,t) .Ui{s,t) 

-ui{s,t)\ T I . ,. .. o^ 



Ul{t,t) 



Ul{t,t) 



Jo- 



(18) 



These identities allow us to remove the initial field configurations 0; , by expressing them in terms of the final fields 
(pf and the derivatives d^± , and obtain the master equation. 

The full equation is very complicated and, as for quantum Brownian motion, it depends on the system-environment 
coupling. In what follows we will compute the diffusion coefficients for the different couplings described in the previous 
Section. As we are solely interested in decoherence, it is sufficient to calculate the correction to the usual unitary 
evolution coming from the noise kernels (imaginary part of the influence action). The result reads 



ipT = (0^|[if,Pr]|<^f ) -iV^T-iDj{uJo,t)pr 



(19) 



where j labels the different couplings. Dj are the diffusion coefficients. The ellipsis denotes other terms coming from 
the time derivative that do not contribute to the diffusive effects. V is understood as the minimal volume inside which 
there are no coherent superpositions of macroscopically distinguishable states for the field. The Tj, which depend on 
the sums and differences ± of the field amplitudes, have to be calculated case by case. 

The effect of the diffusion coefficient on the decoherence process can be seen considering the following approximate 
solution to the master equation 



i"*",^ ;t] exp 



Jo 



(20) 



where is the solution of the unitary part of the master equation (i.e. without environment). The system will 
decohere when the non-diagonal elements of the reduced density matrix are much smaller than the diagonal ones. 

This decoherence time to sets the scale after which we have a classical system- field configuration, and depends 
strongly on the properties of the environment. It is constrained by 



1 « V 



ds Dj{ko,s) > VTi / ds A(fco,s), 



(21) 



for any particular j = I, and corresponds to the time after which we are able to distinguish between two different field 
amplitudes, inside a given volume V. Our conservative choice is that this volume factor is 0{fi~^) since (the 
Compton wavelength) is the smallest scale at which we need to look. 

In the next Section we will compute the corresponding diffusion coefficients, from which we will estimate the 
decoherence time bounds. 



IV. DECOHERENCE TIMES FOR SLOW QUENCHES 

Our earliest results from Ref. [20] were for instantaneous quenches, for which the quench time is effectively /i^^ <C tsp- 
A priori, it was not obvious that this slow quench preserves these results about the decoherence time, that to ^ ^sp- 
This is assumed in the Kibble scenario [5] for domain formation as determined by the appearance of defects, in which 
classical defects are assumed to have appeared as soon as the transition is effected. This was first addressed by us in 
a rather schematic way in [22] , to which what follows provides a more sophisticated analysis and extensions to other 
couplings. 

To tackle this problem, we assume that the quench begins at t = and ends at time t = 2Tq, with Tq 3> ^ f^^^- 
At the qualitative level at which we are working it is sufficient to take m^(To) = fj,^ exactly. Most simply, we consider 
a quench linear in time, with temperature T{t), for which the mass function is of the following form [37] 

{is^ for t < 

^2 _ for < t < 2Tq 

-fM^ " for t > 2rq 

Note that our Tq is the inverse quench rate T~^dT/dt\T=Tct and so differs from that of [37] by a factor of 2. 
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The field behaves as a free field in an inverted parabolic potential for an interval of approximately tgp [36], where 



(22) 



We have rederived tsp, using the exact results of [37] (see Appendix). The end result is that tgp is the solution to 



C 



exp 



4/ Aofap) 
31 t 



3/2' 



(23) 



where Ao(t) =t-Tq^,C = (64v^7r^/^)~\ and t = (rq//i^)^/^ (see details in the Appendix A). 

In order to find the master equation (or more strictly the diffusion coefficient) we follow the same procedure as in 
Section III, assuming a dominant wavenumber fco- The classical solution is given in Eq. (16), where Uj satisfy the 
mode equation 



Ui{s,t) = 0. 



Since we are neglecting the self-interaction term, our conclusions are only believable for t ^tgp. 
The solution of Eq.(24) is given by 



ui{s,t) = 



Ai[-u)^P]Bi[^] - Ai[^]Bi[-u;^P] 



(24) 



U2{s,t) 



Ai[-ulP]Bi[^] - Ai[^]Bi[-ojlP] 



where Ai[s] and Bi[s] arc the Airy functions, with Afe(,s) = ,s — luqP, Ak{t) = t — i^qP, and ojq — fi^ + k^. We note 
that Ao{t) = t — Tq. In the causal analysis of Kibble [5] t (/U~^ <C f <C Tq) is the time at which the adiabatic field 
correlation length collapses at the speed of light, the earliest time in which domains could have formed. Our analysis 
suggests that this earliest time is not t, but tgp. 



A. Quadratic coupling 



We start by considering the biquadratic coupling first, for which the IF is given by 

^2 



ImM = -^Jd^xJ dS A2{x)Nq{x,y)A2{y), 



where Kq{x—y) = ImG'^_|_(a;, y)0{y'^ — x^) is the dissipation kernel and Nq{x—y) = ReG^_^_{x, y) is the noise (diffusion) 
kernel. G++ is the relevant closed-time-path correlator of the field at temperature Tq. We have defined 



1 



/,+ 2 



1 



+ 0-') 



We can formally find the master equation in the same way that we did in Section III. It is given by Eq.(19), with 
a diffusion coefficient of the form (details on how to get this coefficient can be found in Appendix B) 



Jo 



ds u{s, t) F{ko, s, t), 



(25) 



with 



u{s, t) = 



U2{s,t) - 



U2{t,t) 
Ul{t,t) 



ui{s,t) 
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and 



F{ko,s,t) = -i- \ReGl.{0;t- s) + |rcG?_ , (2A;o;t - s) + ^ReG^ , (2^/2A;o; t - s) + ^ReG?. , (2\/3fco; i - s)\ . 

d4 z 4 o J 

It is only in u{s,t) that the slow quench is apparent. G'^_^_{k,t — s) is the Fourier transform of the square of the 
Feynman propagator (x propagator). 

In the high temperature limit (T /i), the explicit expression for the kernels can be shown, with a little labour, 
to be 



ReGl^{k;t - s) = 



Tl_l 
647r2 k 



dp 



C0s[(-\/p2 + ^2 ^ _ g)]^ 



and 



ReG^ (0;t-s) = — ^ / dp , ... 



■cos[2Vp2+/i2(t-s)], 



(26) 



(27) 



where ji is the thermal X-field mass at temperature T ^ Tc- In our scheme, this is approximately the cold x mass. 

It is because the x-^eld propagator is unaffected by the ^-field interactions that the detail of the expressions (26)- 
(27) are possible. Fortimatcly, it is inessential. Wc sec that, for times iit ^ 1, the behaviour of i?qu(fco, t) is dominated 
by the exponential growth of u{s,t), and the integral in Eq.(25) by the interval s ~ 0. Indeed, it is easy to prove 
that u{0,t) = (fflnj)^ ^ 1' "^^il® "(*'^) = 1- Therefore, although the expessions in Eqs.(25), (26), and (27) are 
complicated, we can approximate the whole diffusion coefficient by 



Dq^{ko,t) ^ F{ko,0,t) u{0,t) 



ds 



m(s, t) 
u{0,t) 



(28) 



where we used the fact that F{ko, s, t) is bounded at s = 0. 

We will assume large Afe(i) (and Afc(s)), which means Afe(t), Afe(s) ^ t. This condition is satisfied provided s is 
larger and not to close to a;gTq//x^, and allow us to use the asymptotic expansions of the Airy functions and their 
derivatives for the evaluation of Ui{s, t). This will be justified posthoc. In particular we obtain 

4 fAk{ty 



u{Q,t) ~ 



u(0,f) ~ 



1 



Afc(i) 



1 /Afe(t) 



t 3 



t 



I 4 /Afc(t) 



(29) 



Therefore, it is straightforward to check that in these cases the integral in Eq.(28) is given by 



Wo 



ds 



1 + s 



M(0,i) 
u(0,<) 



Wo 



ds 



and it is 0(1), due to the fact that iit^l. 
We can estimate the decoherence time for the quadratic coupling as 



l>Vr [ dtDqu{ko,t), 
Jo 



(30) 



where V is the decoherence volume, as before. 

In terms of the dimensionless fields (j) = (^^ + (j)~)/2ii, and 5 = (<^+ — (j)~)/2i^ it follows that Pqu = g'^^^cpd'^. 
In order to quantify the decoherence time we have to fix the values of 5, and (p. Inside the volume V we do not 
discriminate between field amplitudes which differ by 0{p), and therefore, for the sake of argument, we take 5 — \. 
For ^ we set = a/A, where A < a < 1 is to be determined self-consistently. It is necessary to be as precise as this 
in the first instance, since tsp and to only differ by logarithms. It is important that the prefactors in the arguments of 
these logarithms are determined carefully, given that terms nominally 0(1) can be small. Once wc have established 
that the difference is large enough, this artificial precision can be dropped. Posthoc it is sufficient to take 6 ~ 0(1) 
and (jp ~ 0(a/A), as we did in [20] and [22]. 
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Then, for quadratic interactions alone the decoherence time reads 

'4 /Afc(tz)) 



exp 



4.10 



4 "^0 



3 V t J [ XT^a' 
Using {(p^) of (A5) we obtain w 7(/i/Tc)(^rq)~^/^. Thus the decoherence time satisfies 



(31) 



exp - 



4 f Akito) 



t 



2.10' 



,3^0 V (MT-q)^ 

2 3 i ' 

/i J.2 ^2 



(32) 



for long-wavelength modes. 

If we compare it with the spinodal time of (23), we get, 



(/iAo(isp))^ - (/iAfe(ir,))i w ^ ^/^Z?^ J 2 In + In 

4 Wo 



1 /7c 

2 U 



(MT-q): 



> 0. 



(33) 



We see that the numerical prefactor in the argument of the logarithm is, indeed, 0(1), and our concern for precision 
was unnecessary, in retrospect. Our approximation scheme depends, as for the instantaneous quench, on the peaking 
of the power in the field fluctuations at long-wavelength fco <C by time tsp, and it is suflicient to take ko « 0. In 
this case (33) becomes 



(MAo(isp))5 - {nAo{tD))^ ~ ^ 



>0, 



(34) 



from which tn < t^p follows, as in the instantaneous case. The inclusion of further interactions, including the self- 
interaction with short-wavelength (fc > /i) modes can only reduce to further. 

We should point out that, whereas peaking is inevitable for the instantaneous quench for weak enough coupling, 
it is not the case for very slow quenches. In such cases our approximations break down and a different analysis is 
required. The details are rather messy, but a sufiicient condition for our approximation to be valid is that yUTq ^ rj/^ 
[36]. Tighter, but less transparent boimds can be given [36]. 

When these bounds are satisfied the minimum wavelength for which the modes decohere by time tgp can be shown 
[22] to be shorter than that which characterises domain size at that time. Although we can talk loosely, but sensibly, 
about a classical domain structure at time tsp wc cannot yet talk about classical defects on their boundaries, as 
the naive picture might suggest. Defects (in this case, walls) are described by shorter wavelength modes {k ^ /x). 
Nonetheless, the classical domain structure is sufficient to determine their density [22]. 



B. Bilinear and linear couplings 

Another possibility that needs to be considered is that of a bilinear coupling of the 0- field to the environment. This 
interaction preserves the cj) ^ —cj) symmetry of S'syst['^] in which, for simplicity we continue to take bilinear couplings 
equal, as g'^ = g' / \/N. We treat this interactions as an additional set of interactions to the biquadratic interactions, 
whereby is qualitatively unchanged. 

The influence action is still obtained from Eq.(12), as 



ReJ^biiin = — 7— 



J d^x J d% A2{x)Kb{x - 2/)E2(y), (35) 



and 



ImMbilin = / '^^^ / "^^y ^2(a;)iVb(a; - y)A2{y), (36) 

where 

/^b = ReG++(a;,y)%0-a;0), N^, = lmG++{x,y). (37) 
The temporal diffusion coefficient is now given by 
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-Dbiiin(fco,t) = / ds u{s,t) Fuiin{ko,s,t), (38) 
Jo 



where 



Fbiiin{ko,s,t) = |lmG'++(0;t - s) + h.mG++{2ko;t - s) + hmG++{2V2ko;t - s) + hmG++{2V3ko;t - s)| . 



with 



ImG++(fc;s) = -^—^^—JLA. (39) 



Following the same arguments as in the previous section, we can evaluate the diffusion coefficient by 



DbMko, t) « Fbiii„(fco, 0, t)^^M, (40) 

Wo 



where u{Q,t) is given in Eq.(29). Thus, 

£'biiin(/co, t) « 10-2^^ u(0, t). (41) 
A* Wo 

For the unstable long-wavelengths {k^ < approximately) we find, for times > 1, that -Dbiiin(fco5 *) again 
shows the exponential growth 



I5bm„(fco,i)«^^V^exp|-(^^j |. (42) 

It follows that Fbiiin = As a result, the contribution to the decoherence time from bilinear interaction (for 

long-wavelength modes) can be evaluated as 

The value of a is again determined from the condition that, at time to, {(j)^) ~ arf. That is, = 0.1(/LtTq)~^/^. 
Thus, decoherence time is given by 



From (23), to and t^^ are related by 

3 3 3 r 1 

(/iAo(tsp))^ - {ii^o{tD))'i « ^ ^/M^ In [3(MTq)4J , (45) 

which is positive for a sufficiently slow quench. 

As we have already observed, early studies of decoherence were confined largely to quantum mechanical systems, 
for which the environment was typically a collection of harmonic oscillators, to which the system coupled linearly. 
Such systems have the virtue of exact solvability (or closed equations) and have been very instructive. However, in 
the context of quantum field theory linear terms are usually a signal of an inappropriate choice of field basis. Further, 
as we are considering models with spontaneous symmetry breaking, linear couplings with external fields are not a 
natural choice since they break the vacuum degeneracy. We include isolated linear couplings for completeness. Again, 
choosing couplings equal, as g'l = g" /^fN, and defining Ai = (0+ — (/>~)/2 and Si = (0+ + (/>~)/2, we are able to 
write the real and imaginary parts of the infiuence functional as 

ReMiin = J d^x J d'^y A^{x)Kbix - y)^i{y), (46) 

and 
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(47) 



For times nt ^ 1 the diffusion contribution to the master equation is 

Ain(fco,*) « \/w(0,t). 

Note that the exponent is only half that of the quadratic and bihnear interactions. 
In this case, it can be shown that the decoherence time comes from 



(48) 




5 5 _2 l 



(49) 



where r[a, z] is the incomplete Gamma function and, in this case Fun = \g"^fi^d'^. As to '> t, we can approximate 
Eq.(49) and obtain, 



exp • 



j~ 9 11^ 11^ 



(50) 



The decoherence time associated to the linear interaction term, can be written as 



Afe(tz3)5«^P <!31n— +ln^^^ 
^ ^ ' M 9 yti* Tc ^5 



16 r?^ 1 oj§ n'^ 



(51) 



to and tsp are related by 



(/iAo(isp))^ - {fJ,Ao{tD))^ ~ - Va*^ In 



10^^(/xr,) 



(52) 



In this case, due to the bound on the quench time (rq < r]/j2^), we find that it looks as if tzj > tgp- However this 
result is not believable as it stands, since the diffusion coefficient has been computed assuming t < tsp. Whatever, the 
rapid decoherence of the biquadratic and other couplings is not present. This shows how adopting linear coupling to 
an environment in mimicry of quantum mechanics can be misleading. 

C. Comparison between different couplings 

In our present model the environment fields Xa are not the only decohering agents. The environment is also 
constituted by the short-wavelength modes of the self-interacting field 0. Therefore, we split the field as = +</>>, 
and define the system (</><) by those modes with wavelengths longer than the critical value A^^, while the bath or 
environment-field ((^>) contains wavelengths shorter than A~^. In order to consider only the unstable modes inside 
our system, we will set this critical scale A of the order of /z. In practice, where the separation is made exactly is 
immaterial [36] by time tjj, when the power of the 0- field fluctuations is peaked at fco ^ The effect is to give a 
separation of system from environment through the decomposition of S[(j),x] of (2) introducing an interaction terms 



of the form 



4>t(t>t, and 



t>%4'>- The relevant interaction term is 



5'couple[</><,X] = -Sint [(/><, X] - ^ j d'^x{(t)<{x)(j)y{x)f. 



(53) 



All terms omitted in the expansion [28,29] do not contribute to the one-loop calculations for the long-wavelength 
modes that we shall now consider. ^ 



^Strictly speaking, for yu/3 < 
See Ref. [28] for details. 



< yU one should include an additional term proportional to ( 



in the interaction Lagrangian. 
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The net consequence of the separation of the long-wavclcngth system modes of (j) from the short-wavelength 
modes (/)> of the environment through the interaction action ^couple of (53) gives an additional one-loop contribution 
D^iko, t) to the diffusion function with the same u{s, t) and the same form as in (25). However, G±± is now constructed 
from the short-wavelength modes </>> of the (/)-field as it evolves from the top of the potential hill. Without the 
additional powers of A''"^ to order contributions the one-loop calculation is unreliable. In fact, we would not expect 
the inclusion of the (^field to give a qualitative change at one-loop level. In the first instance, the approximation 
(28) remains valid, and the 1-loop term is driven by the exponential growth of u{Q,t). Secondly, the effect is that 
the short- wavelength modes have been kept at the initial temperature Tq, on the grounds that passing through the 
transition quickly has no effect on them. That is, with .g ~ A and no factor, the short- wavelength (f) modes give a 
contribution comparable, qualitatively, to all the explicit environmental fields put together. At an order of magnitude 
level there is no change, since the effect is to replace g'^ by + O(A^) = 0{g'^). 

However, since the contribution of D^{kQ,t) to the overall diffusion function is positive we can derive an upper 
bound on the decoherence time to from the reliable diffusion functions _Dqu(fco,t), £'biiin(A:o, i) (and I?iin(fcoj i))- Let 
us suppose we have a theory where the three considered couplings with the environment are present. Eq.(21) will be 
satisfied because one of its terms will grow faster than the others, rather than because many terms will each give a 
small fraction of unity. Specifically, we have (up to numerical factors 0(1)) 

/•to rto rto 

Tqu / ds £'qu(fco, S) : rbilin / ds £)bilin(fco, s) : / ds Diin{ko,s) 

~.:j(t):4(t)-{im'}- 

Since n and 5 « 0(1) then, for jito ^ 1, we have 

f^D r^D r^L 



r^D r^n r^D 

Tqu / dsZ)qu(fco,s) > rbilin / rfs £>bilin(fco, «) > / ds Ain(fco, «), 

Jo Jo Jo 



(55) 



if g,g' and g" are comparable. Thus, it is sufficient to evaluate the constraint on the time to from this biquadratic 
interaction alone. 



V. FINAL REMARKS 



We first summarize the results contained in this paper, in which we have shown how the environment leads to the 
decoherence of the order parameter after a transition. After the integration over the scalar environment-fields Xa in 
Section II, we have obtained the coarse-grained effective action (CGEA) for the system (field (p). From the imaginary 
part of the CGEA we obtained the diffusion coefficient of the master equation, at 1-loop and in the high temperature 
environment limit. Terms omitted are relatively 0(-/V~^/^) for N weakly coupled environmental fields. Subsequently, 
we evaluated the decoherence time for the long- wavelength modes of the system-field (for each of the different couplings 
with the environment considered) for slow quenches. This decoherence time depends on the coupling between system 
and bath, the self-coupling of the system (through the environment temperature T), and the mass. In our model, we 
have shown that the decoherence time is in general smaller than the spinodal time. 

We stress that the inequality tu < tgp is insensitive to the strength of the couplings, for weak coupling. The reason 
is twofold. Firstly, there is the effect that F cx Tq , and Tq oc is non-perturbatively large for a phase transition. 
Secondly, because of the non- linear coupling to the environment, obligatory for quantum field theory, T rx (jr. The 
completion of the transition finds (jr ~ rj"^ oc A~^ also non-perturbatively large. This suggests that F, and hence to, 
can be independent of A. In fact, the situation is a little more complicated, but the inequality holds. 

This arises because the diffusion coefficients, which trigger classical behaviour when they become large enough, are 
controlled by the exponential growth of the unstable modes. It is this same exponential growth that determines tgp. 
This result provides a post-hoc justification of our assertion [36,38] that the spinodal time sets the scale for the onset 
of classical behaviour (in open systems). 

Our emphasis has been on the many weak environments because of the control that this gives us on establishing a 
robust upper bound on to- However, we noted earlier that their total contribution was qualitatively that of the short- 
wavelength modes of the (p field alone. Environmental fields are an important feature of the early universe, but even 
had we not included them we would have expected a similar result from the one- loop couplings of short-wavelength to 
long-wavelength (p modes. Although, in this case, we have no way to control the higher loop terms, it is quite probable 
that our prescription of the environments is unnecessarily detailed, and early decoherence is a general feature. 
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Ideally, we could extend our ideas to gauge theories. This is beyond the scope of this paper, but we would like 
to emphasize an important point. As diffusion is additive (each extra term in the interaction action gives additional 
diffusion coefficients to the master equation), the inclusion of further couplings and fields can only reduce to. Further, 
the coupling of short to long-wavelength modes of the global theory is omnipresent. However, theories with derivative 
couplings (scalar QED, for example) tend to produce small additional diffusive effects at low k (scaled by k). It is for 
this reason that, if scalar fields are present, there is no need include gauge field interactions in order to estimate an 
upper bound to the decoherence time. 

The results obtained in this paper also justify in part the use of phenomenological stochastic equations to describe 
the dynamical evolution of the system field, as we will now discuss. As it is well known [28,29], for (jy^x"^ interactions 
one can regard the imaginary part of SA as coming from a noise source ^mix), with a Gaussian functional probability 
distribution. Taking the biquadratic coupling to the external environment x^Seld first leads to a noise, termed ^2, 
say, with distribution 



P[6(x)] = TVj, exp - 



1 



(56) 



where N^^ is a normalization factor. This enables us to write the imaginary part of the influence action as a functional 
integral over the Gaussian field ^2{x) 



V^2{x)P[^2] exp 



i{A2ix)Ux) 



exp { - I I dS: I d^y A2{x) g^N^{x,y) A2{y) 



In consequence, the CGEA can be rewritten as 



where 



A[<t>+,ct>-] = -Un j 2?6i^[6]exp|zSeff[</.+ ,r,6]|, 



5eff[</.+ ,0-,6] = ReA[cP+,<j,-] - J d^x[A2{x)^2{x) 



Therefore, taking the functional variation 



0, 



we are able to obtain the "semiclassical-Langevin" equation for the system- field [28,29] 

a(l){x)-fi^(f>+^cf>''{x)+g^cj>{x) j d^yKq{x-y) c^^y) = ^{x)^2{x), 



(57) 
(58) 
(59) 

(60) 

(61) 



where jl and A are constants "renormalised" because of the coupling with the environment. Since (61) is, from (60), 
the stationary phase approximation, it is only valid once the system has become classical. It can be used to establish 
domain formation only because to < tsp- 

Each part of the environment that we include leads to a further 'dissipative' term on the left hand side of (61) with 
a countervailing noise term on the right hand side. Once we include the interactions between ^< and at one-loop 
level we get a similar equation to (61) from the term S'coupio alone (sec Eq.(53)). However, although the 0^ip> 
4>'l^4>y terms were ignorable in the bounding of in the Langevin equations they give further terms, with quadratic 
(/)^^3 noise and linear noise ^1 respectively. 

The inclusion of bilinear interactions leads to the inclusion of further terms in (61) of the same form. For the linear 
interaction with the environment (to the exclusion of self-interaction) we do recover the additive noise that has been 
the basis for stochastic equations in relativistic field theory that confirm the scaling behaviour of Kibble's and Zurek's 
analysis [25]. 

For times later than tsp, neither perturbation theory nor more general non-Gaussian methods are valid. It is difficult 
to imagine an ah initio derivation of the dissipative and noise terms from the full quantum field theory. In this sense, 
a reasonable alternative is to analyze phenomenological stochastic equations numerically and check the robustness of 
the predictions against different choices of the dissipative kernels and of the type of noise. We stress that this is only 
possible because to < t^p. This will be considered further elsewhere [39]. 

Finally, we see that the role initially attributed by Kibble [5,23] (and subsequently by others e.g. [40]) to the 
Ginzburg regime is just not present. 
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APPENDIX A: THE SPINODAL TIME 



In this Section, we show the estimation of the spinodal time tgp defined from {(p'^)t=tBp ^ rf 
The equation of motion for the mode Uk{t), with wavenumber k is, in the quench period. 



1 

H s 



Uk{t) = 0, 



(Al) 



subject to the boundary condition Uk{t) = e for t <Q, where 



+ fc2. 



Instead of the simple exponentials of the instantaneous quench, Uk (t) has solution 

l(k{t) = akAi[^]+bkBi[^], 



(A2) 



with Ai[s], Bi[s] the Airy functions; Afc(f) = t — lo^P and t = (rq//u^)^/^. Note that Ao{t) = t — tq, the time since 
the onset of the transition. In the causal analysis of Kibble [5] i <^Tq) is the time at which the adiabatic 

field correlation length collapses at the speed of light, the earliest time in which domains could have formed. Our 
analysis suggests that this earliest time is not t, but tgp. 
It is straightforward to establish a relationship between t and tsp > i. The constants of integration in (A2) are 



ak = wlBii-oj^f) + iojtBi{-uj^f)], 



(A3) 



It follows that, when Ak{t)/t is large, then 



■ cot 



jit 



Ak{t) 



Ao(i) 



1/2 



exp 



1/2 



exp 



t 



3/2-1 



3/2-1 



(A4) 



where P = t-^/^{Ao{t))-^/^ /2. 

For large initial temperature Tq = 0{Tc), we find the power spectrum for field fluctuations peaked around k, and 



2w 



[ k^ 



dk \Uk{t)[' 



CTofAojt) 



-5/4 



exp 



AoW 



3/2- 



(A5) 



We have intentionally included the prefactor C to show that terms, nominally 0(1), can in fact be large or small (in 
this case C = (64-\/27r'^/^)~^ = O(10~^)). Note that, although the unstable modes have a Hmited range of fc-values, 
increasing in time, this is efi^ectively no restriction when Ao{t)/i is significantly larger than unity. 
Finally, we obtain 



Tc 

— ~ ^ exp 



ij 
C 



4 Mo (tsp) 
3 V t 



3/2- 



(A6) 



where C = C[ln{np7]'^ /CT^)-^/^]. Since the effect on tsp only arises at the level of 'In In' terms, C" « C is a good 
estimation in all that follows. Since this choice underestimates tsp it only strengthens our results that tgp > to- 
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APPENDIX B: THE DIFFUSION COEFFICIENT 



In this Appendix we show how to obtain the diffussion coefficient Eq.(25) (for the quadratic couphng) from the 

reduced propagator. Following the same techniques used for the quantum Brownian motion to obtain the master 
equation we must compute the time derivative of the propagator Jr (Eq.(14)), and eliminate the dependence on 
the initial field configurations (j)f' that enters through the classical solutions (jt^^ This can be easily done using the 
propagator identities of Eq.(18). Thus, using Eq.(16) and identities (18), wc can write the initial field configuration 
4\ , in terms of the final configurations {(l>f), final field derivatives {d^±), and the functions of time Ui{s,t), 

id^t - cl>.) m{s,t) Jo = u,{s,t) - cl>i) Jo + (Bl) 

neglecting terms proportional to derivatives respect to the final field configuration which do not contribute to normal 
diffusion. 

The temporal derivative is given by 

lhdtJr[<Pt^'t^iA't>t^^iA = |/tren['/'+]-/lren[r]-»g' ^'^^ ds A^'(s) F{ko;S,t) 

where i^q(A:o; s, t) is the Fourier transform of the dissipation kernel. The ellipsis denotes other terms coming from the 
time derivative which do not contribute to diS^usion. 

Difi^usive effects are associated with terms proportional to — in the master equation. Using Eq.(Bl) we 

remove initial conditions from Eq.(B2), and looking only those terms proportional to , wc get the master equation 
Eq.(19). The diffusion coefficient D{ko,t) comes from the noise contribution to the influence functional and it is the 
time dependent coefficient that multiplies A^^. Thus, for the quadratic coupling example, we find 

Dqu{ko,t)= [ ds u{s,t) F{ko;s,t), (B3) 
Jo 

where 



u{s,t) 



U2{s,t) 



Ul{t,t) 



ui{s,t) 



For different couplings between system and environment one should follow the same procedure shown in this 
Appendix. 
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